BD-A148  711 
UNCLASSIFIED 


MULTIDIMENSIONAL  INVERSE  SCATTERING  FOR  FIRST  ORDER  1/1 

SVSTEMS(U)  CLARKSON  UNIV  POTSDAM  NV  INST  FOR  NONLINEAR 
STUDIES  A  I  NACHMAN  ET  AL  JUL  84  AFOSR-TR-84-0918 
AF0SR-78-3674  F/G  1271  NL 


AD- A 148  711 


TR-  84.0918 


INSTITUTE 

FOR 


I.F.N.S.  #38 


f 


NONLINEAR 

STUDIES 


MULTIDIMENSIONAL  INVERSE  SCATTERING  . 

FOR  FIRST  ORDER  SYSTEMS 

Approved  for  I 

^  distribution  uali*il84* 


Adrian  I.  Nachman*  and  Mark  J.  Ablowltz 
Clarkson  University 


Clarkson  College  of  Technology 
Potsdam,  New  York  13676 

♦Permanent  Address:  Department  of 
Mathematics,  University  of 
Rochester,  Rochester,  NY  14627 


DTiC 


DEC  0  6 1964 


July  1984 


r *■»  »“  ri 


84  10  28  30 


■*»  »*■  »V»lii**.  «V»*viV«Vt 1  i*fi\iS.i*i,il»l‘i,,s 


K>-AMriL 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


rt  CLASSIFICATION  AUTHORITY 


OSSIFICATION /  DOWNGRADING  SCHEDULE 


JRMING  ORGANIZATION  REPORT  NUMBER(S) 


AME  OF  PERFORMING  ORGANIZATION 

Clarkson  University 


6b  OFFICE  SYMBOL 
(If  applicable) 


ADDRESS  (Gfy.  State.  and  ZIP  Coda) 

Dept  of  Mathematics  &  Computer  Science 
X  Potsdam  NY  13676 


I  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

)  AFOSR 


i  ADDRESS  (City,  State,  and  ZIP  Code) 

!  Bolling  AFB  DC  20332 


8b  OFFICE  SYMBOL 
Of  applicable) 

NM 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 

unlimited. 


S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

AFOSR- TR-  8  4.0 


Ja.  NAME  OF  MONITORING  ORGANIZATION 

Air  Force  Office  of  Scientific  Research 


7b.  ADDRESS  (C/ty,  State,  and  ZIP  Code) 

Directorate  of  Mathematical  &  Information 
Sciences,  AFOSR,  Bolling  AFB  DC  20332 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

AF0SR-78-3674 


TO  SOURCE  OF  FUNDING  NUMBERS 


*53 

IV 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO 

NO. 

NO 

ACCESSION  NO 

61102F 

2304 

A4 

title  (Include  Security  Classification) 

MULTIDIMENSIONAL  INVERSE  SCATTERING  FOR  FIRST  ORDER  SYSTEMS 


-^;z  PERSONAL  AUTHORS) 

Adrian  I,  Nachman*  and  Mark  J.  Ablowitz 


I 


1 3a.  TYPE  OF  REPORT 

Technical 


13b.  TIME  COVERED 
FROM _ TO 


14.  DATE  OF  REPORT  (Year.  Month.  Day) 

JUL  84 


>  j6  SUPPLEMENTARY  NOTATION 

V  *Permanent  address:  Dept  of  Mathematics,  Univ  of  Rochester,  Rochester  NY  14627. 


It.  SUBJECT  TERMS  (Continue  on  reverse  If  necessary  and  identify  by  block  number) 


-  ■* 

COSATl  COOES 

FIELO  i 

I  GROUP  1 

1  SUB-GROUP 

9.  ABSTRACT  (Continue  on  reverse  If  necessary  and  identify  by  block  number) 

A  method  for  solving  the  inverse  problem  for  a  class  of  multidimensional  first  order  systems 
is  given.  The  analysis  yields  equations  which  the  scattering  data  must  satisfy;  these 
equations  are  natural  candidates  for  characterizing  admissible  scattering  data.  The  results 
are  used  to  solve  the  multidimensional  N-wave  resonant  interaction  equations. 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 


»>>  Qunclassifico/unumiteo  □  SAME  AS  RPT  □  OTIC  USERS 
NAME  OF  RESPONSIBLE  INDIVIDUAL 
s  Dr.  Robert  N.  Buchal 


ftv  LJ  FORM  1473, 84  MAR  *3  APR  action  may  ba  used  ui 

■  •  All  Afhtr  mem  n 


21  ABSTRACT  SECURITY  CLASSIFICATION 


202)  767- 


83  APR  adition  may  ba  utad  until  aihauftad 
All  othar  aditioni  are  obwlata. 


22b  TELE°HONE  (Include  Area  Code)  J  22c  OFFICE  SYMBOL 

NM 


m 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 

- 


KrtflAfrgcr  vvwMfirjAMkiiSik  ammsitf  -n-co&oco 


i 


V  .  I  " .  *  _  ■  .  •.  «.  *1 


Accession  For 

MULTIDIMENSIONAL  INVERSE  SCATTERING  FOR  FIRST  ORDER  SYSTEMS'  rn.;  cR'ii 

■’  i  r:  c  T'  •: 


Adrian  I.  Nachman*  and  Mark  J.  Ablowitz 
Department  of  Mathematics  and  Computer  Science 
Clarkson  University,  Potsdam,  New  York  13676 


ABSTRACT 


A  method  for  solving  the  inverse  problem  for  a  class  of  multidimensional 
first  order  systems  is  given.  The  analysis  yields  equations  which  the  scatter¬ 
ing  data  must  satisfy;  these  equations  are  natural  candidates  for  characterizing 
admissible  scattering  data.  The  results  are  used  to  solve  the  multidimensional 
N-wave  resonant  interaction  equations. 


1.  Introduction 

The  inverse  scattering  problems  for  the  hyperbolic  and  elliptic  generaliza¬ 
tions  in  the  plane  of  the  m  x  m  AKNS  system  have  been  successfully  studied  in  [1] 
and  applied  to  the  linearization  of  several  physically  significant  nonlinear 
evolution  equations  (N-wave  interaction,  Davey-Stewartson,  etc.)  in  two  spatial 
and  one  temporal  dimensions. 

We  indicate  here  how  the  method  used  in  our  investigation  of  the  n-dimensional 
Schrodinger  equation  [2]  can  be  applied  to  the  study  of  the  inverse  problem  for 
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La  f  3xq  +a  JJi  Sx^  '  Q(x0’x^  *  ^ 

Here  are  constant  real  diagonal  m  x  m  matrices  (we  denote  the  diagonal 
entries  of  by  and  assume  t  J]|  t  0  whenever  i  t  j);  the  matrix¬ 
valued  off-diagonal  potential  Q  =  (Q1J)  may  depend  on  xQ  as  well  as  x  =  (x-j . xn) 

and  a  =  +  iaj  is  a  complex  parameter. 

The  operator  (1)  is  associated  with  the  nonlinear  system: 

lOli  =  1  a  lOll+Wa  ji  B*\  9Qlj  +  1  _a  io\ 

3t  o  aij  3xq  J'a1jJ£  V  3x^  a  2(aU  a£j'Q  Q 
(with  a. .  =  -4 - f  ,  1  <  l  <  n,  for  some  real  1  <  l  <  n,  1  <  i  <  m). 

•J  _  I1  ~  *  ~ 


Even  though  no  traditional  scattering  operator  exists  in  the  case  <jj  +  0,  the 
so-called  3  method  (see  [2]  and  references  given  there)  gives  a  satisfactory 
definition  of  scattering  data  for  L^,  along  with  a  systematic  inversion  pro¬ 
cedure,  which  we  use  to  solve  (2). 

A  solution  of  the  inverse  scattering  problem  for  the  hyperbolic  case 
Oj  =  0  is  then  obtained  by  a  limiting  argument;  thus  we  can  avoid  a  separate 
study  of  a  Riemann-Hilbert  boundary  value  problem  (which  in  higher  dimensions 
may  also  involve  some  geometric  complications).  The  main  advantage  of  this 
approach  is  that  it  yields  (from  the  compatibility  conditions  associated  with 
3  in  several  variables)  equations  which  must  be  satisfied  by  the  scattering  data. 

In  addition  to  their  importance  for  the  problem  of  characterizing  admissible 
scattering  data,  these  equations  have  several  consequences:  i)  they  provide  a 
formula  for  reconstructing  the  potential  from  the  scattering  transform  purely 
by  quadratures  (in  the  generic  case  when  no  three  of  the  vectors  J1  *  (jj  .J^,. . .  ,J^), 
1  <  i  <  m,  are  col  inear);  ii)  they  show  how  one  can  recover  the  scattering 
transform  from  (at  least  small)  data  given  on  certain  (n+1 )-dimensional  surfaces 


(n+1  being  the  number  of  variables  in  Q);  iii)  they  may  indicate  what  other 


(possibly  non-local)  evolution  equations  could  be  solvable  with  the  1ST 
developed  here;  iv)  they  constitute  in  themselves  a  new  class  of  multidimensional 
nonlinear  systems  of  integro-differential  equations  which  are  linearizable. 

2.  The  Case  t  0. 

We  will  denote  by  k  =  (k^,...,kn)  =  kR  +  ikj  a  point  in  Cn  and  will  often 
write  f(k)  instead  of  f(kR,kj)  for  an  arbitrary  function  of  kR  and  kj. 

As  a  first  step  in  the  §  procedure  we  construct  a  family  of  solutions  of 
L  ^  =  0  of  the  form  ip  -  y(xQ,x,k)exp[ii£_i  ^(x^-oxqJ^)]  with  y  bounded;  y  will 
then  satisfy  the  equation 

5^ +  oZ"=i  Jn  lx^ +  k*cVu] ' Qu-  (4) 

The  generalized  eigenfunctions  y^  »  (y^J)  we  will  work  with,  are  obtained  by 
solving  the  integral  equation  ua  »  I  +  G0(Qya),  i.e. 


uoj  =  6ij  +  fj  Gaj(xo'yO’X"y,k^Q^yO’y)ua(yO*y’k))ljdyOdy’ 


where  the  Green's  function  is  given  by 


i(xft£n  +  x*0 


For  brevity  we  will  assume  here  that  Q  is  such  that  this  integral  equation  has 
a  bounded  solution  ya  for  all  k  e  Cn. 

Ga  can  be  computed  explicitly  by  contour  integration: 


4i  slgn(o.jJ)  ia^(x0,x,k)  n 

V  (xn»x»k) - - e  °  °  n  6( 

0  0  ZTri^-ojjxg)  1-2 


\ 4  *  A  df  “  o  X« 

ao  (xo,x,k)  s  ^  a  ( 1°  I  Xrtkr  -  —j"  (ak^),). 
a  o  ls]  ax  0  x,  1 

The  next  step  is  to  express  Sy  in  terms  of  y.  We  start  by  writing 


with 


(8) 


-4- 


a 

M 


i 


rCi 


'vp  JSp 

Z—  and  hence  -=—  (Qu)  as  a  superposition  of  exponentials: 

3k  3k 

P  P 

3G  o(j!  -  jj)  r  n  .  ig13(xn,x 

M-  (Qu  )),J  -  - 4  °<  1  JlVe  ° 

3k0  0  2i  |aj  |  (27r)n  ^n  4=1  11 


n  •  iB  3(xn,x,k,X) 

o(  I  Jlxje  a  0 


T^(k,X)dX.  (9) 


.  .  «  .  II  • 

with  e^J(x0,x,k,x)  =  a^J(xn,x,k)  +  2  (Xa-CpjJxJX,,  and 


■  u  R  £A0'  X, 

A”  1 


'oj(M>  .  ff 


"iBaJ(y0’y,k,X* 


( Q (yQ .y )ua (y0 .y » k ) )  J  dyQdy 


The  calculation  of  3y  is  then  based  on  the  following  crucial  symmetry 
property  of  our  Green's  function: 

“i^fTJ(xn’X»k»X)  r.:  ri  ^i  i  i 

e  U  G^J(x0,x,k)  =  G^(x0,x,k^J(k,X))  whenever  SJ^X^O;  (12) 

here  k^(k,X)  is  the  point  in  Cn  whose  £th  component  is 

jJ  -  j1 

(kaJ(k»X))jj,  =  ~~1  -  (aki)l  +  ki  +  Xsl  *  (13) 

Once  (12)  has  been  established  it  can  be  shown  (assuming  that  (5)  admits 
no  homogeneous  solutions)  that 


y  o(J-JJ)  r 

=2-  (xQ,x,k)  -  2  - - -  6(2J>  )T^J(k,X)i 

k  0  i  ,j  2i  |oj  | (2tt )n  ^n  1  0 


i®QJ (xQ’X’k’X) 


*  Ma(x0»x,k^(k,X))EijdX; 


(we  have  denoted  by  the  m  x  m  matrix  with  entries  E^  =6^ r<S js ) •  If  we  now 
fix  all  k^,  i  f  p,  and  apply  the  (1 -dimensional )  inhomogeneous  Cauchy  integral 
formula 


to  the  kp  variable,  we  can  convert  (14)p  to  an  integral  equation:  noting  that 
u(xQ,x,k)  -v  I  when  |kp|->-»  (and  denoting  k*  =  (k.| ,. . .  ,kp, . . .  ,kn) )  we  have 


ua(*0,x,k)  =  I  - 


i  i  fff  i i  i (Xg,x,k 1  ,A) 

i!d(Jp-Jp  )Jll^rT-(kU)e 


x  yo(x0,x,Kij(k,,X))E..dAdkJ  dkj  .  (16) 

J  p  p 

(More  generally,  one  can  use  (15)  with  f(z)  =  Pa(xg»x,k+zv) ,  zed,  with  k  fixed 
and  with  an  arbitrary  v  e  dn  which  is  not  perpendicular  to  any  of  the  vectors 
i  f  j).  The  matrix- valued  function  Ta(k,A)  defined  in  (11)  is  our 
scattering  data  and  (16)  is  the  inverse  scattering  recipe  for  reconstructing 
y  from  T.  Once  y  is  found,  the  potential  is  easily  recovered: 


Q(x0.x)  =  10  [Jp,  j|  _0  (xQ,x,k)dkR  dk,  ]. 


r  3*p  P  P 

On  the  other  hand,  given  an  arbitrary  T(k,x),  to  apply  the  above  inversion 

procedure  we  would  first  need  to  know  that  the  equations  (14)  ,  p  =  l,2,...,n, 

2.2  p 

are  compatible;  requiring  that  v  yields  the  following  characterize 


3k  3k  3k  3k 
r  p  P  r 


tion  equations  for  T: 
•i  -i  _  4  3T 


Lpr[V  *  +  4  (Jp-Jp)(J)-jjr)(jV^-jr^)  = 


1*1  i 1 

x  T^V^k.X'),  X  -  Vx')dx'.  .  08: 

Q  i  ■ 


For  compatibility,  (18)^  need  only  hold  whenever  £j!a„  ■  0,  however  one  may  also 


verify  that  Ta  when  given  by  (11)  satisfies  (18)  everywhere. 

It  turns  out  to  be  very  useful  to  recast  (18)  in  integral  form.  It  is 

enough  to  keep  only  the  equations  (18)pl •  We  then  look  for  a  parametrization 

of  the  hyperplane  ((k,A)  e  Cn  x  Rn  :  =  0}  by  new  variables  (x»Wq,w)  e  <jn~^  x 

x  R  x  Rn  so  that,  in  the  new  coordinates  L1}  =  -4—,  2  <  p  <  n  and 

p‘  3Xp  =  = 

6^J(xo»x.k,X)  =  XgWg  +  x*w;  these  requirements  determine  (up  to  a  translation 
of  x)  the  following  (invertible)  map: 


k£  ”  (Jl  ~  Jl^xt*  1  =  2;  kl  ~  +  j1  i  wo+Z!=lJtwJl^; 

J^-J1  |a| 

(jj-0])(jj-jj)  4 

■  77  +  v*i 2>  xr£^[-  .V--  (B**>i  - 


Vi 


Vi 


(19) 


•  • 

To  use  (15)  as  before,  we  need  the  limit  of  T1J  for  |y  |  large  (and  y.,  *  t  p, 

P  * 

Wg,  w  fixed);  this  turns  out  to  depend  on  whether  for  some  r  f  i , j  we  have 


(20) 


For  brevity  we  consider  only  two  cases  (the  only  ones  arising  in  the  study  of 
(2)  -  see  the  appendix):  case  I-relation  (20)  does  not  hold  for  any  distinct  i,j,r 
and  any  p  f  1 ;  and  case  II-relation  (20)  holds  for  all  i,j,r,p  (in  other  words,  the 
vectors  J Jm  all  lie  on  the  same  line  in  Rn).  In  the  generic  case  I  we 
have 

TaJ(x»wo’w)  =  Q1j(w0,w) 

IXpl-*00 

and  08)p]  becomes 


(21) 


-7- 


)1J(w0,w)  =  l 
■ 


-1(x0w0  +  x-w)  ,J 


Q  (*q»x)  dXgdx  and  \  (X2* • • • *Xp» • • • »X^) 


If  (20)  holds  for  some  r  f  j  then  (21)  need  not  be  true  (see  (7),  (8),  (11)). 

gyl  j 

In  case  II  we  have  -  =  0  for  all  p,  2  <  p  <  n;  this,  together  with  Li ouvi lie's 

3Xn 


theorem,  allows  us  to  replace  (22) j  by 
TaJ(X‘Vw)  =  T’j(wn*w) 


a  v"0’ 


In  case  I  we  conjecture  (as  in  [2])  that  the  main  condition  needed  to 
characterize  the  scattering  data  is  that  Ip^[To](x»w0,w)  be  independent  of 
X  and  p  and  have  suitable  decay  properties  in  (Wq,w);  furthermore,  given  a  T^ 
which  passes  this  admissibility  test  we  can  (re)construct  a  local  potential  Q 
simply  as  the  inverse  Fourier  transform  of  I [T] . 

From  (22) n  we  see  that  T1J  is  completely  determined  by  its  values  on  the 
(n+1 )-dimensional  surface  x  =  Xq;  the  analogue  of  this  in  case  I  is  the  follow¬ 
ing:  given  T^J(xo»wo»w)  =  G1j(Wq,w),  1  %  i,j  <  m  we  have  (from  (22)j) 

ii  ii  iff  *Vj[T  ](x'  »wn»w)  ^Dl  ^a^Xn»w0’w) 

t>,w0,w,  ■  .’>vo  ♦  -  pV07--?p-°-]Vy-ir 


p  p 


which  (at  least  for  small  G)  could  be  solved  to  find  T  everywhere. 


3.  The  case  a  =  -1 . 

If  we  formally  substitute  a  =  -1  in  (6)  we  find  that,  away  from  the 
hyperplanes  E.j  =  (k  e  tn  :  =  the  eigenfunction  u_-j(xQ,x,k) 

is  well-defined  and  holomorphic.  Thus  it  appears  that  the  inverse  problem  for 
the  hyperbolic  system  L_-|  could  be  regarded  as  a  Riemann-Hilbert  problem  with 
data  on  the  hyperplanes  1  <  i  <  j  5  m.  We  prefer  to  obtain  an  inversion 

procedure  from  our  results  for  Oj  t  0.  There  seems  to  be  little  advantage  in 


studying  the  limit  of  ya(xQ*x»k)  as  a  -1  (it  leads  us  back  to  an  analysis 
of  singularities  on  the  hyperplanes  Z.  .);  we  work  instead  with  the  limit  of 

*  J 

ua(Xg,x,kR,ajkj) ,  with  kj  now  playing  the  role  of  a  parameter.  From  (6)  we  find 
(Xg  »x»kR  >Ojk  j )  —  (xg ,  x » kR » k  j ) 


(2iOn+1 


I'CXqCo  +  X*C) 

X  0  u  u 


d£0d£. 


(24) 


with  e( • )  the  Heaviside  function;  correspondingly,  lim  y  (xg,x,kR,Ojkj)  = 

a  -1+iO 

=  yL(xg,x,kR,kI)  where  yL  solves  the  integral  equation  yL  =  I  +  GL(QyL).  From 
(24)  we  see  that  yL(xg,x,kR,kj )  is  a  solution  of 


^  '  z"=l  Jl  $7  '  1£*-1  ■  Qli 


(25) 


for  every  value  of  the  parameter  kj  in  F  .  Our  scattering  data  is  now 

(xn*x’^R»ki  »A)  •• 

L  U  X  i  /n/„  /„  „  L  L,  \  \  1 J , 


TLJ(kR,krx)  = 


n+1 


(Q(^0»x^L^x0,x,kR’kI^  dXgdx 


(26) 


.  .  .  .  X 

with  (xg>x,kR,k j ,X)  *  zj=1 JjfJji^xOkI  '  T^R  "kI  ^  +  ^x£+Jix0^Aji^  Takin9 

1  J£  1  1 

limits  in  (14)  we  find  the  reconstruction  equations  for  y: 


yL(x0’x,kR’kI )  =  1  +  (2u)n+1  iEj (JP_Jp\  J.  ^kR  +10  +  Fr  -k‘  -?0  ^6(iJjiV 


0 (kj  -kj  )  e(kj  -kj  ) 


ii  i$|(x«»x,kR,kj,X)  ^ 

(kR’kI’A)e  yL(x0,x,^J(k',kj,x))Ei.dxdk^  dkj  , 

J  r  p 


X 


(27 


where  now  (k  J(kR,kT,x))R  s4k.  +  -  •  ■  kT  +  X.  and  (k.  J  )T  =  kT  . 

L  K  1  K£  y  j1  l  L  li 

&  9+ 

To  write  the  characterization  equations  for  we  introduce  new  variables 
(suggested  by  the  limit  of  (19))  (xr»Xi»wq»w)  e  to  parametrize  the 

hyperplane  Zjjx^  =  0  in  F3n  as  follows: 


kR  ■  (Jj  -  J])XR  .  A  >  2;  kR  =  eJ=2(J^  -  JJ)XR  +  -pj  (VZ5=1J£w£] 

£  £  1  £  J  ~  J  -j 

I  =  2;  ki  =  zJ=2(J£"4)xi  +  TT3  wo  ( 

£  £  I  £ 


(VXi^r  A  >  2;  X1  =  Z”a2[ . ! . . . - 


(V\] 


Then  under  the  assumptions  of  case  I  in  section  2,  the  limit  of  the  equations 
(22) j  yields: 

0(xi  -xj  )  e(xi  -X!  )  .. 

TLJ(Xr’Xi>Vw)  =  Q1J(w0,w)  +  ^||CXr"xr  ho  +  XR  -Xr  -i0  Kl[TL](x' *w0’w)dxR  d*I 


P  P 


P  P 


with  Np-|[J|_]  given  by  a  slight  modification  of  (18).  In  case  II  we  have 

^Xr»Xj>wo»w)  =  (wq * w ) *  (29)j 

As  in  section  2,  we  can  now  use  (29)j  to  characterize  admissible  TL,  (recon¬ 
struct  Q,  as  well  as  recover  T^  from  data  given  on  xr  =  const.,  xj  =  const. 

It  should  be  pointed  out  that  once  the  family  of  Green's  functions  has 
been  chosen,  all  the  above  results  can  be  derived  without  recourse  to  our  limit¬ 
ing  arguments  (vk  uL  can  be  expressed  in  terms  of  uL  using  the  appropriate 
symmetry  property  of  G^  and  the  analytic  behaviour  of  for  kj  large  -  needed  to 
establish  (27)-  follows  from  (24);  these  analytic  properties  together  with  the 


yyyyy.v. % v -av.a^v , 
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compatibility  requirements  ^ ^ —  imply  (29)) 

^p  ^p  V 


4.  Relation  between  T.  and  the  Scattering  Operator  (o=  -  1 


To  fix  notation  we  sketch  an  elementary  definition  of  the  scattering 
operator  associated  with  L  -j .  When  Q  =  0,  given  f  :  Fn  +  Rm,  the  solution  of 
the  Cauchy  problem  L_^u(xq,x)  =  0,  u(0,x)  =  f(x)  is 

u^Xq.x)  =  f1  (x-j+XqJ^  , . . .  ,xn+xQj|!)) ,  1  <  i  <  m,  which  we  write  as  u(xQ,x)=f  (x+xQJ) 

When  Q  is,  say,  smooth  and  of  compact  support,  given  any  (reasonable) 
f  :  Fn  -*■  Fm  there  is  a  unique  u  solution  of  L_^u  =  0  with  u(xQ,x)  =  f(x+xQJ) 
for  xQ  «  0;  furthermore  there  is  a  unique  g  such  that  u(xQ,x)  =  g(x+x^J)  when 
Xq  »  0.  We  write  g  =  Sf.  On  the  Fourier  transform  side  S  can  be  written  as 

Sf(S)  =  f(C)  +  — —  S(?,kR)f(kR)dkR.  (30) 

(2ir)  JRn 

The  question  we  would  like  to  address  is  how  to  recover  T^  (and  hence  Q)  given 
S(C,kR).  To  relate  T^  and  S(£,kR)  it  turns  out  that  we  need  to  relate  uL  and 
the  eigenfunction  y(xQ,x,kR)  corresponding  to  the  "Volterra"  Green's  function 

Glj(x0,x,kR)  =  Q(xg)exp[-iE^_i  (x£+xo^l^R  ^  n  5^x2,+xOJ^'  (31) 

We  start  with  the  identity 

UL  -  y  =  (Gl-G)(Qul)  +  G(Q(yL-u)),  (32) 

write  G^J-G1J  as  a  superposition  of  exp(isj^)  and  use  a  suitable  symmetry  pro¬ 
perty  of  G.  The  main  result  is  the  following  linear  equation  for  TL  given  S: 

T?J(kR,kT,X)  =  Sij(kjj(kR,k  ,X),kR)  -  — 2  f  0(2  jj\)  x 
L  R  1  R  R  I  R  (27T)n  £=1  X.  Z 

Sii'(kjj(kR,kI^)rR,j(kR,kI,X'))T^:i(kR,kI,Xl)dX',  (33) 


'  »  *  .  -  *■  •  .  *  •  *  .  •  .  •  •  •  »•  ■  *  •  *  I  ■  *  «  ’  »  •  •  .  *»’.*.  ’  *  *  -V*  «.*  V*  *  •  •  *  .  '  ■  (L*  «  ' 


x 
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where  k^J(k^,kj,A)  stands  for  the  real  part  of  k^ 


ij 


5.  Applications  to  Nonlinear  Equations 

The  equations  (2)  are  the  compatibility  conditions  (cf.  [  ])  for  the 

Lax  pair: 

Lo*  '  0  and  ft  <1  B*  ^  <*» 

the  matrices  B^,  1  <  i  <  n,  are  constant  real  diagonal  and 

A1'’(t,XQ,x)  =  1  a^ (t,XQ,x)  with  a.j  given  by  (3).  The  restrictions  imposed 

by  (3)  on  the  matrices  J  ,  1  ^  %  <  n,  are  discussed  in  the  appendix.  To  find 

the  time  dependence  of  the  scattering  data  corresponding  to  (2),  we  set 

^  =  y  exp[i  ^ ^ x£_axo^ fi.) ^ *  then  ^  sat"'sfles  (^)  as  we^  as 


Au  #  It  +  z"=i  B£  ISt  +  iz"=i  k*[BrM]  ■  Ay  =  °- 


(35) 


Applying  the  operator  A  to  both  sides  of  the  equation  (14)  we  find  (when  Oj  t  0) 


9T 


ij 


3^“  (t.k.X)  =  i2j=1[BjkrB^J(k,A)]T^(t,k,A). 


(36) 


For  the  case  a  =  -1  the  equations  (obtained  as  limits  of  (36)  or  by  a  parallel 
derivation)  are 

3T.lj  .  ... 

-ft-  (t,kR,kj»X)  =  l£j=1[BjkR  -BjkJJ(kR,kI,X)]TLJ(t,kRfkI,X).  (37) 

Thus,  when  a  =  -1 ,  we  can  apply  the  inverse  scattering  procedure  together 
with  (37)  to  construct  the  solution  to  the  Cauchy  problem  for  (2).  Note  that 
TL(t,kR,kj ,A)  as  given  by  (37)  satisfies  the  characterization  equations  if 
TL(0,kR,kj ,X)  does. 

When  Oj  ?  0  the  Cauchy  problem  for  (2)  is  ill-posed;  however  (by  analogy 
to  the  corresponding  linear  problem)  we  can  use  inverse  scattering  to  solve  (2) 
as  follows:  given  Q(0,xq,x)  it  can  be  decomposed  into  Q+(0,Xg,x)  +  Q_(0,Xq,x) 


^  A 


WVV 


where  Q+{0,xq,x)  extends  to  a  function  Q+(t,Xg,x)  satisfying  (2)  in  the  half¬ 
space  t  >  0,  while  Q_(0,Xq,x)  extends  to  a  function  satisfying  (2)  in  the 
half-space  t  <  0.  Assume  for  simplicity  that  aTa.  .  >  0  for  all  i  f  j.  Given  Q 

1  I  J 

A  —  A 

define  Q+  by  Q+(0,Wg,w)  =  0(+Wq)Q(O,Wq,w) .  If  T+  is  the  scattering  transform 
of  Q+  then  from  the  direct  problem  we  find  (0,x»Wq,w)  =  0  for  wQ  >  0;  thus 
for  t  >  0  we  can  define  (see  (36))  (t.x.Wg.w)  by 


Tj.j(t,x,w0,w)  =  exp[itzJ=1(Bjk£-B^k^)]T|j(0,x.w0,w)  =  (see  (3),  (13)  and  (19)) 

=  exp[it(-iiwQ  +  zJ=1(a1j.jJ-B^)w£)]T|j(0,x,w0,w).  (38) 


Since  the  expression  in  the  exponential  does  not  depend  on  x  and  since  its  real 


part  is  nonpositive  if  t  >  0,  (t.x.Wg.w)  satisfies  the  character!' zati 


on 


equations  (29)  so  inverse  scattering  should  yield  the  desired  potential 
Q+(t,Xg,x) ;  similarly  we  can  construct  Q_(t,Xg,x)  solution  of  (2)  for  t  <  0. 


Acknowledgements 

It  is  a  pleasure  to  thank  A.S.  Fokas  for  helpful  conversations.  This 
work  was  supported  in  part  by  the  National  Science  Foundation  under  grant  number 
DMS-8202117,  the  Office  of  Naval  Research  under  grant  number  N00014-76-C-0867, 
the  Air  Force  Office  of  Scientific  Research  under  grant  number  AF0SR-84-0005  and 
by  a  University  of  Rochester  faculty  research  grant. 


jqj 

cCcC 


Appendi x 

We  need  to  find  trie  restrictions  imposed  on  the  choice  of  matrices 
1  ^  l  S  n,  by  the  existence  of  (a^j)  and  Bf,  1  s  £  S  n  satisfying  (3). 

If  (3)  holds  then  the  matrix  (a.  .)  is  symmetric  and 

jj-jj 

3ip  "  apj  =  ^ai j-apj ^  ~p~T  for  Gvery  1  and  every  i»j»P  distinct  .  (A1 

J£_J£ 

(Conversely,  if  (Al)  holds  with  (a-jj)  symmetric  then  B^,  1  <  £  <  n  can  be 
found  so  that  (3)  is  satisfied.)  Mote  that  if  i  a  y  (Al)  forces  J1 ,  JJ, 

JP  to  be  col  inear.  There  are  two  cases: 

I  a.  =  a  .  for  all  i,j,p  distinct.  Then  (Al )  puts  no  restriction  on  J  ; 

I P  PJ  J6 

in  particular  they  can  be  chosen  so  that  (20)  does  not  hold  for  any  distinct 
i,j,r  and  p  f  1 .  The  system  (2)  is  linear  in  this  case. 

II  For  some  i n , jn , pn  distinct  a-  f  a  .  .  We  show  that  in  this  case  the 

1  m  0P0  P°J° 

vectors  J  , . . .  ,J  must  all  be  col  inear.  From  (Al )  we  already  know  that 

po 

J  ,  J  ,  J  are  col  inear.  For  any  r  f  io’^0’P0  0ne  the  ^°"'"'ow’in9  mu?,t  be 


(i)  a.  f  a  •  ,  (ii)  a.  ^  a .  „ 

V  rJo  r,o  Vo 


(i i  i )  a  .  f  a  . 

rj0  JQP0 


(for  if  not  a.  =  a  .  =  a  .  =  a  •  contradicting  our  assumption).  In  either 

'opo  no  rjo  poJo  f0  p0  j- 

of  the  possibilities  (A2)  J  will  be  on  the  line  passing  through  J  ,  J  ,  J  ; 

this  will  be  true  for  any  r,  1  <  r  <  m.  (Conversely,  given  j\j^,...,jni 

col i near  with  f  we  can  construct  ( j )  symmetric  satisfying  II  and  (Al)). 

12  m 

It  follows  that  whenever  (2)  is  not  linear,  the  matrix  having  J  ,J 

as  rows  has  rank  at  most  2;  thus  if  n  >  3  its  columns  (the  diagonals  of  the  matrices 

J  in  (1))  must  be  linearly  dependent  and  then  the  inverse  scattering  problem  for  L 
X/  o 

can  also  be  solved  by  reducing  it  to  a  lower  dimensional  one.  On  the  other  hand, 
since  the  characterization  equations  are  trivial  (i.e.  N(T)  =  0)  in  this  case, 
it  seems  reasonable  to  expect  that  other  (possibly  non-local)  nonlinear  equations 


can  be  found  which  would  be  compatible  with  (22)^ 
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